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PREFACE 


Two separate phases of propellant sloshing under low-gravity conditions were investigated during this 
project. Results of the first phase, an experimental and analytical study of sloshing in rectangular tanks, 
have been reported previously in Technical Report No. 1, January 1970. A summary of Technical Report 
No. 1 is presented as the Appendix to this Final Report. Details of the second phase, which is comprised of 
an analytical study of sloshing in ellipsoidal tanks, are included as the major portion of this Final Report. 

The SwRI Project Manager was Dr. Franklin T. Dodge. Mr. Luis R, Garza contributed substantially to 
the experimental portion of the research, and Dr. Wen-Hwa Chu made several valuable suggestions con- 
cerning the analytical approach. 

The entire program was made possible by the continuing efforts of the NASA-MSFC technical 
monitors: Mr. Robert S. Ryan, Mr. Frank Bugg, and Mr. Harry Buchanan. 
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ABSTRACT 


A variational integral is formulated from Hamilton’s principle and is proved to be equivalent to the 
usual differential equations of low-gravity sloshing in ellipsoidal tanks. It is shown that for a zero-degree 
contact angle the contact line boundary condition corresponds to the “stuck” condition, a result that is due 
to the linearization of the equations and the ambiguity in the definition of the wave height at the wall. 

The variational integral is solved by a Rayleigh -Ritz technique. Results for slosh frequency when the 
free surface is not “bent-over” compare well with previous numerical solutions. When the free surface is 
“bent over,” however, the results for slosh frequency are considerably larger than those predicted by 
previous finite-difference, numerical approaches; the difference may be caused by the use of a zero degree 
contact angle in the present theory in contrast to the nonzero contact angle used in the numerical 
approaches. 
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I. INTRODUCTION 


The study of liquid sloshing under strong capillary and weak gravity conditions, usually called 
low-gravity sloshing, was begun in detail with the publication in 1964 of the report [1] * by Satterlee and 
Reynolds, dealing with circular cylindrical tanks. By 1966, low -gravity sloshing in a circular cylindrical tank 
was sufficiently well understood to allow an equivalent (mathematical) mechanical model to be 
derived [2] . Detailed research concerning spherical and ellipsoidal tanks began somewhat later, with the 
delay probably caused by the knowledge that there would be considerable analytical difficulties with these 
geometries. Thus, it was not until early in 1968 that the present analytical work was started at SwRI to 
predict low-gravity sloshing behavior in an ellipsoidal tank, under the initial sponsorship of Contract 
NAS8-20290, and later under the sponsorship of Contract NAS8-24022, both with NASA-MSFC. 

To avoid as much as possible the expected analytical trouble in satisfying exactly all the matliematical 
“boundary conditions” for low-gravity sloshing in spheroidal tanks, it was decided to formulate the equa- 
tions of motion as an integral of the kinetic and potential energy (Hamilton’s principle) and then to find 
the solution by a Rayleigh-Ritz approximation technique. Rattaya had shown this procedure to be satis- 
factory for ordinary sloshing (flat interface) in spheroidal tanks [3] . In fact, Yeh had already derived such 
an integral for conditions under which the free surface does not “bend over” although he presented no 
numerical solutions [4] . 

During the 3 years since this research was started, several numerical solutions of the exact differential 
equations of low-gravity sloshing have appeared, either for an ellipsoidal tank [5] or for a more general 
axisymmetric tank [6] .f Thus, the numerical results of the present work are now somewhat out-of-date, 
especially since a completely satisfactory solution for sloshing with a “bent over” free surface has not been 
obtained. The present analysis is still significant, nevertheless, because several difficulties are pointed out 
which relate to the difficulty in obtaining correct expressions for the behavior of the slosh wave at the tank 
wall in the event that the contact angle is zero degrees. These difficulties are not so evident in a completely 
numerical treatment, where the contact angle is usually taken to be about five degrees [5, 6]. Further- 
more, an approximate but easily made calculation of the equilibrium free-surface shape is presented herein; 
this result should be useful for making quick computations of liquid volumes or tank sizing requirements. 


* Numbers in brackets refer to References listed in Section VI of this report. 

f A major reason for the delay in completing the present work was the Project Manager’s heavy involvement in developing 
the computer program described in Ref. 6. 


II. EQUILIBRIUM FREE SURFACE 


When gravitational forces are small, the shape of the free surface is determined by the combined 
influence of gravity g, surface tension a, and contact angle \p. Thus, the equilibrium free surface for an 

axisymmetric tank is a surface of revolution about the 
z-axis (the z-axis and the direction of gravity coincide) as 
shown in Figure 1 . 



FIGURE 1- NOMENCLATURE FOR EQUILIBRIUM 
SHAPE OF LIQUID 


Let the equation of the free surface be z=/(r) and 
the equation of the tank shape be z = w(?-); both /and w 
may be double -valued functions. The free-surface shape 
is determined by the following equations [1 ,2] : 


- p = aC 


( 1 ) 


which relates the pressure jump to the curvature C of the 
free surface. The curvature C is 


C 




ridfldr) 


/■ dr l[l + (d//dr)2] 


1/2 j" 


( 2 ) 


where the negative sign corresponds to the bent-over 
part of the surface. The liquid pressure p must satisfy 
the hydrostatic requirement: 


P=Po-PSf 

where p'^ is the liquid pressure at r = z = 0. Combining Eqs. (1), (2), and (3) gives 


( 3 ) 


^ a d r r(d//dr) 

"7di^l [1 + (d//dr)2]i/2 


- Pgf=Po-Po 


( 4 ) 


The pressure difference pg —p'g depends on the curvature at the center ( r - z - 0). Because/ - d//dr - 0 at 
r = 0, Eq. (4) shows that this pressure Jump is 




r^O 


( 5 ) 


Thus, the governing differential equation can be obtained by combining Eqs. (4) and (5): 


a d 
+ — - 


ridfldr) 


rdr \[l + idfldrf]^^^ 


■ 2a 


dr^ 


-pg/=o 


r = 0 


( 6 ) 


The boundary conditions on / are 


df 

f= — = 0 at ?• = 0 
dr 


( 7 ) 



4/ + arctan 



atr = rw 


( 8 ) 
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where ;// is the angle of contact. By letting R = F = f/ro, W= and Eqs. (6), (7), 

and (8) can be written in nondimensional form as 


d r R{dFl6R) 

) /d^F' 

R dF \[\ +(dF/dRy] 

F=0 1 


dF 
dF ” 

atF = 0 

1 

ii 

atF =Fyv 


-NboP'=0 

R = 0 


( 9 ) 

( 10 ) 

( 11 ) 


and, in Eq. (11), the contact angle \jj is assumed to be zero. The tank shape is h^(/^) = —H ±L{\ 
where H - h/r^^ and L = (see Fig. 1). 

The nonlinear Eq. (9) cannot be solved in closed form. Furthermore, there were no computer 
programs to solve Eq. (9) numerically for bent-over free surfaces at the time the present work was started, 
although several now are available. Thus, an approximate solution was obtained. After a considerable 
amount of trial-and-error analysis, the following equations were found to represent the free-surface shape 
adequately. For the lower part of the surface, 


^’ = Fi (12) 

and, for the upper, or bent-over, part oi the surface, 

F = =/35 [1 -(7?/5)"] (13) 

The point at which Fj is joined to F 2 R- B- b/ro , which is the nondimensional radius to the point 
where the surface has its vertical tangent (see Fig. 1). For surfaces that are not bent over,F 2 is not needed, 
and, in this case, B may be larger than Fw The parameters )3, ju, and B are determined by the boundary 
conditions on F and by geometric considerations. 


If the surface is not bent over, then Eq. (1 1) is satisfied if 


|3 = TF 




^ 1/2 


1 -R 


2 

W 


Requiring that the free surface must intersect the tank wall atF = gives the result that 


(14) 


[\-{\-{RjBfVI^^=-H-L{\-Riyi^ (15) 

Finally, Eq. (12) must satisfy Eq. (9) at least on the average; that is, F, must satisfy Eq. (9) when 
integrated from F = 0 to F = F. This requires that 


1 

& 


= 1 +0.168F^A^bo 


(16) 


Equations (14), (15), and (16) allow j8, B, and Fy^ to be calculated for a given Nqq and volume 
of liquid. 
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If the surface is bent over, the same procedure shows that 


nRl, 



_ 1 -Rl 


( 17 ) 


+ =-H + L{\ 


(18) 


1 +0.1685^A^bo 


(19) 


Furthermore, the curvature [Eq. (2)] of both Fj and F 2 must be the same at F = F in order to give a 
continuous pressure. Thus, 


( 20 ) 

2 

These equations specify Fj and F 2 for a bent-over free surface. 

Comparisons of the free-surface shapes predicted by these equations to the exact shape given in 
Ref. 5 for a 5° contact angle, are shown in Figures 2a, 2b, and 2c for a spherical tank and various Bond 
numbers. The comparison is very good except for the case of a nearly full tank. 
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AXIAL COORDINATE. 




a. Nbo = 10. 75% and 12.5% Full 


0- Ngo = 30, 87.5% and 25% Full 



e. Nbo = 100, 62.5% and 37.5% Full 


FIGURE 2. COMPARISON OF APPROXIMATE AND EXACT THEORIES 
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III. SLOSHING ANALYSIS 


A. Basic Equations 

The integral or variational form of the sloshing equations of motion can be obtained either from 
Hamilton’s principle or by a direct integration of the differential equations of motion. In both cases, it is 
necessary that the “kinematic” boundary condition for the wave height be satisfied independently; that is, 
77 and 0 must satisfy the relation that 

077 _ 30 d/ /30\ _ 30 
bt dz dr \3r / dn 

where 77 , the wave height, is measured vertically above the equilibrium free surface and d(p/dn is the fluid 
velocity normal to the free surface. This equation is nondimensionalized by assuming that 0 = 0[r,0,z] 
cos u)t - cos and 77 = 7 ? [7^,0] sin cjt = (rQ/Q)e [R, 9] sin r2r, where = cj(r^/^)^/^ is 

the nondimensional natural frequency and r = nondimensional time. In nondimensional 

form, the above equation is thus 

e = 3^/3Z”(di^/dR)(3$/3/^) = (34)/3/7)[l +(dF/d/^)^] , Z = F(R) (21) 

According to Hamilton’s principle, the remaining equations of motion can be obtained by requiring 
that 0 and 7 ? are such as to minimize the integral 

I = j (T-U)dt 


1 + 


dr 


1 - 


, z=Ar) 


where T is the kinetic energy of the liquid and U is the potential energy, including surface tension effects, 
and the integral is taken over one cycle of the motion. Following S a tterlee and Reynolds [1] , this integral, 
in nondimensional variables and after the integration over time is performed, is 

27T 


/ = f fj (y^-V^)RdRde dZ-^(^j^'^J 
V ^ 

4 iMA 


(deldRy 

[1 +(dFldR) 


[1 +{dF/dRf]^l'^ 

(deIRdey 


do 


\R=R, 


[1 +{dFldR) 


2i 1/2 


}••') 


RdRdO 


( 22 ) 


where the first integral is taken over the entire volume V of the liquid, the second integral over the contact 
line R = Ry^^, and the third integral over the free surface^. The authors of Ref. 1 apparently believed that 
Eq, (22) does not lead to the correct differential equations because, after formulating an integral similar to 
Eq. (22), they then use a different variational integral in their numerical work. However, Eq. (22) does in 
fact yield the correct equations, as will now be shown. By proceeding as shown by Hildebrand [7] , the 
variation of Eq. (22) is 


2tt 


6/=n" J J ■ 8V^) dV -(R^riNBo) J 


ede 


[1 +{dFldRy]^l^ 


dd 


\R=R, 


. , / 1 f ideldR)8{deldR) ^ jdeldd) 8 (deldO) 1 , \ 

l[l +(dF/di?)2]3/2 R^[\+(dFldRy]'l^ J ^ / 


RdRdO 


(23) 
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Because 6(V^) = the volume integral can be transformed by Green’s tlieorem to give 


Vd^)dV=- JJJ (V^^)8<^dV+ ^‘** ‘^ ”**^ 

V V y 


where the first surface integral is taken over the free surface?^ and the second over the tank walls V. The 
surface integrals in Eq. (23) are now integrated once by parts to give the following relations [7] : 



(de/dR) 8 (de/dR) 
[1 +(dF/dRf]^l'^ 


^ 2tt 

I i? di? d0 = J ■ 


R (deldR) 6e 

[1 +(dF/dR)^p/^ 


R =R^ 


R = 0 


dO 


T 


TJ_. 

r R (deldR) '1 

1 RdR 

L[1 +(dF/d7?)^]3/2 J 


8eR dR dO 


(24a) 


and 


//{ 

T 


(de/dd) 5 (de/dd) 


R^[l +(dF/dRy] 


2il/2 




-//{ 

7 


RdRdd 

0 

R^[\ +(dP/d/?)2]^/2 


(ae/30)6e 


[1 +(dF/d/?)2] 


RdRdd 


2il/2 


0 = 27T 


R dR 


6 = 0 


(24b) 


The first integral on the right in Eq. (24b) is equal to zero because 5e at ^ = 0 is the same as 5e at 0 = 2ti; 
that is, 8e is continuous. Collecting terms, 


6 / = - 


///<’= 


4>)64>dK + ^2^ 


II 


— 1 S^dT^^- 

dn 




II 


1 


R (de/dR) 


+ ■ 


it)l 

0 

d'^e/dO^ 


Nqo l[l +(dF/dR)^]^^^ j ' R'^[\ -\-(dF/dR)^]^^^ 

1/2 


[(3e/a7? ) - re] 5e 
[1 +(dF/d/^)^]^/^ 


de 


R^R^ 




a4> 

dn 


1 + 


dR 


54> 


RdRdd 


(25) 


where the fact that the element of area is [1 + (dF/dR)^ ] ^ i? dR d0 has been used in the last term. 

Because 64> and 6e are completely arbitrary, 81 can only be equal to zero (that is, the value of the 
integral can be minimized) if each of the integrals in Eq. (25) is separately equal to zero: 


= 0 in the liquid volume, V 
d^ 

— =0 on the tank walls, Z = W(R) 
dn 


and 


3e 

dR 


Te = 0 on the contact line ,R = R^ 


(26) 

(27) 

(28) 
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Equations (26), (27), and (28) are consequences of requiring that the first three integrals on the 
riglit-hand side of Eq. (25) be equal to zero. The last integral in Eq. (25) involves both 5$ and Se; since the 
relation between 5<h and 6e on the free surface is not known, setting the integrand in this term equal to 
zero does not lead to a useful result. But, if in the original variation, the volume integral is transformed in 
an alternative way as 


Iff 




=-/// 


(6V^4>)4>dF + 


II 


3^ Y 

dn ) 


$dF+ jy 

■w 


a<f>Y 

dn ) 


^dlV (29) 


then a useful form does result. According to Eq. (26), the variation can be zero only if = 0; thus, we 
also have 6V^<h = 0. Further, according to Eq. (27), d^/dn = 0 on”W ; therefore, 6(d4>/3«) = 0 onlM . And 
because d^/dn = e[\ + (dF/df?)^ ] “ ' onTandF dF d0 = [1 + (dF/dF)^ ] - ^ d^^, the surface integral 
over T on the right-side of Eq. (29) reduces to// <^8eR dR d6. Using this last relation and then proceeding 

T 

with the rest of the manipulations as before, the variation of the integral in this alternative form is 


5/ = -f2^ 


Iff [6V«4.] 

V 




[(de/dR) - Fe] 6e 
[1 +(dF/dF)^]3/2 


■dd 

R =i?vv 


rrfJL 

3 

r R (dejdR) 

1 d^e/de^ 

JJ \Nbo 

1 

1[1 +(dF/dF)2]3/2 , 

[1 +(dF/dF)^]‘/^ 




-e + 


deRdRdd 


(30) 


Consequently, the only way that 5/ can equal zero is if in addition to requiring that Eqs. (26), (27), and 
(28) be satisfied, it is also required that 


1 

r a j 

r R (de/dR) 1 

[ d^e/d9^ 

^BO 

RdR 1 

L[1 +(dF/dR)^]^>^ J 

I [1 +(dF/dF)^]*/^ _ 


, Z-F(R) (31) 


Eq. (31) is the last of the required sloshing equations of motion. Thus, any function <I>, with e calculated 
from by Eq. (21), that minimizes the integral Eq. (22) will also satisfy all of the sloshing equations of 
motion. 


As mentioned previously, Satterlee and Reynolds [1] and Yeh [4] used a variational integral that was 
derived by a direct integration of Eqs. (26), (27), (28), and (31), with Eq. (21) as a required side condition. 
This procedure leads to 

27T 2 




Nbo I 


(de/dRy 


- + 


(deldOy 


[\ +(dFldRy]^l^ R^ [1 +{dFldRy]^l^ 




+ eH i? dR d0 



(e^)R dR de 
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(32) 


Eq. (32) differs from Eq. (22) in the sign of the first surface integral and in the addition of the term 
^ d0. There is no evident physical significance to Eq, (32), but it does lead directly to the 
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differential equations of motion, which can be proved by repeating the procedure that yielded Eq. (25). In 
this case, all of the Eqs. (26), (27), (28), and (31) fall out at once, without further arguments or repeating 
the variation in an alternate form. For this reason, Eq. (32) may lead to a more rapid convergence of the 
sequence of trial functions for to the true solution; thus, Eq. (32) or a similar formulation will be used in 
the present work. 


Although Eq. (32) is suitable for use with a free sur- 
face that is not double-valued or bent over, the linearization 
of the equations with respect to e breaks down at the point 
on a bent-over free surface where dF/dR = ^ (i.e., at the 
point of vertical slope of the surface). At this point, e = 
according to Eq. (21), unless also d^/dR ~ 0. But 3^/dR = 
0 is not a realistic requirement unless the point where the 
slope is vertical is also the point where the free surface 
meets the wall; in this case d^/bR ~ d<^/dn does equal zero 
because the normal velocity at the wall must be zero. In 
general, however, e = and thus the second and higher 
powers of e cannot be neglected in the equations. On the 
other hand, this prediction of an infinitely large wave height 
is simply a consequence of the way in which the wave 
height is defined. If instead of measuring e vertically above 
F, e is measured in the direction normal to F, as shown in 
Fig. 3, then Eq. (21) is replaced by 

e = d^/dn at Z = F(R) (33) 



FIGURE 3. NORMAL-TANGENTIAL 
COORDINATE SYSTEM 


For this case, e is always finite. (As will be seen presently, however, this definition of e still leads to 
difficulties in the contact angle condition.) To be consistent, Eq. (32) also must be modified to account for 
the change in the definition of e: 


27T 


I = fff (V^-V^)dV-n^ ff (e^)dr-^(R^r’lNso) f e 

(de/duy 


dO 


R =R, 


r f ^ f / de \ 2 1 

jJ U +(d^’/d^)" \Rde ) ^ j 


[1 +(dF/dR) 


2n 1/2 


dT (34) 


where 


^ (d^F/dRy , (dF/dRf 

(j — I T I ■ 


[1 +(dFldRyy R^ [1 FidFldRf] 


(35) 


The double-valuedness in the analytical representations of F, and e along the free surface can be removed 
by transforming the coordinate F,0,Z = F(R) which locates a given point on the surface into the coordinate 
S,0, where S is the arc -length along the surface from the point /? = Z = 0 to the point in question.. The 
transformation is 


— = [1 +(dF/cy?)2]>/2 
dR 


Z = F(R) 


(36) 


as shown in Fig. 3. 
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The normal derivatives, 3d>/d«, on the free surface and the tank >valls can be evaluated by the 
relations 


3^ 

dn 

bn 


”3(F 


_3Z 

W j 


3^ 

,/dH^\ 

3Z 



)-]/ 


[1 +(dF/d^)^] 

[1 +(dW/dRf]^^^ 


Z = F(R) 
, Z = 


(37a) 

(37b) 


where the change in sign is required because the positive normal direction must always point outward. 


Finally, for numerical work, it is best to eliminate the volume integral in Eq. (34) by requiring that all 
of the assumed functions for must satisfy = 0 identically. Then, Green’s theorem shows that 
Eq. (34) reduces to 




+ 


II 




3<J> 


dV--l 


2 

\dn / 2 

J J 

T 


{deimf 


JJ \Nbo 

7 

_1 + {dFIdRf ^ 

\Rde, 


2 \Nbo J 


-Ge^ 


+ - 


[1 +{dJ^ldRY 


R =Rv 


1/2 j" 


dd 




(38) 


An analogous equation can be written when e is measured vertically above F. 


B. Contact Line Condition 


(22) or (32) is 
Fe 2 LtR=R^ 


r'e 2itR~Ry^ 

For the “no stick” condition (the liquid slides freely along the wall) of interest here, F and F' are functions 
only of the tank geometry and the free-surface shape, and the magnitude of F or F' is picked so that the 
contact angle that the wave makes with the wall always is equal to the static contact angle. But, F or F' 
cannot be calculated unambiguously because it is impossible to define a wave height e at the wall in an 
unambiguous manner. This is demonstrated in Fig. 4, which shows the point of the wave contacting the 
wall cannot be correlated in a definitive way with any point or; the static-free surface, either when e is 
measured vertically above F or when measured normal to it. As can be seen, various values of F or F' are 
possible, depending upon what kind of approximations are used. The approximation used by Concus, et 
al. [5] and by Chu [6] for the normal direction definition of e (Fig. 7c) leads to 

1 r d^WjdR^ cos (d"F/d7?") ) 

sini// 1 [ 1 +(db//dF) 2 ] 3 /^ [ 1 +(dF/dF) 2 ] 3/2 J ’ 

and thus, F' = when \jj = 0. This requires that e = 0 at for otherwise be/bS - which implies 

that there is a jump in the wave height at the contact line. 


The contact line condition implied by Eq. 

bR 

and the one implied by Eq. (34) or (38) is 

bS^ 
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Now,, b^l^n = 0 all along the walls; 
and, for the case v// = 0, on the 

wall at = R\v 3lso equals on the 

free surface at R = /?vv - Therefore, b<i^/bn 
on the free surface at the contact line also 
is zero; but, because e = b^jbn on the 
contact line (and all over the free sur- 
face), it can be concluded that e = 0 on 
the contact line, which agrees with the 
first conclusion concerning e at 
In fact, e ~ ^ ?lX R = is the “correct” 
contact condition when \jj = 0 and the 
equation be/bS = F'e is essentially redun- 
dant. But, the fact that e = 0at/?=/^vv 
requires analytically that the liquid is not 
allowed to move away from the equili- 
brium contact line; in other words, the 
linearized slosh wave is stuck'' to the 
contact line. This is probably not true in 
practice but is just a mathematical con- 
sequence of both the linearization of the 
equations and of the ambiguity in defin- 
ing a proper wave height at the wall. 


EXTENSION OF FREE SURFACE 



Ref. 5 


In the finite-difference numerical ^t \7 

analyses presented in Refs. 5 and 6, the 

contact angle \p was always equal to a FIGURE 4. VARIOUS CONTACT LINE CONDITIONS 

non-zero value (usually 5 degrees). Hence, 

in these studies, P' and e 0 at the contact point. (Further, b^jbn onT at the contact point is not 
equal to b^/bn onTM at the contact point, and thus e does not have to be zero at = Rw) H might be 
possible, therefore, to approach the case \p ^0 numerically in order to determine if the numerical solutions 
converge to a definite value as \p ^ 0. In this regard, however, Chu [6] noticed a sensitivity in his results to 
changes in the value of \p or P^ for small values of \}j. \n the absence of evidence to the contrary, it seems 
likely that the definition of P^ in Eq. (39) is not sufficient to eliminate the mathematical difficulties when 
= 0, and, for this reason, it might not be accurate to extrapolate the results presented in Refs. 5 and 6 to 
the case \p = 0 unless the extrapolations are verified by accurate experiments. 


Measuring e vertically above F also leads, in general, to the result that e = 0 at the contact line when 
\p = 0. The only exception to the condition that e = 0at/^=i?^is the case for which the free surface at the 
contact line has a vertical tangent, as in a cylindrical tank; for this case, the vertical wave height can be 
defined unambiguously at the contact line. (Mathematically, dF/dR at R -R^y ^t\d so e = (3<F/9«)[1 + 
(dF/dRy] = 0 X which is undefined. Thus, e is not necessarily zero at R = R^, and, in fact, the 
correct contact line condition is bejbR = 0 at R = R^.) 
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IV. SOLUTION OF EQUATIONS 


A. Rayleigh-Ritz Technique 

The form of Rayleigh-Ritz technique used here to solve the variational integral involves assuming a set 
of functions = 1 ,2,3 . . each of which satisfies = 0. Then, the velocity potential is 

written as the sum 

m,0,Z)= Y O‘n’^r,(R,0,Z) 

n = 1 

where the constants are to be determined. The wave height e, calculated by Eq. (21) or Eq. (27), is 

N 

e(R,d)= ^ a„e„(R,d) 
n = 1 

These expressions are substituted into the variational integral, Eq. (38), and the various definite integrals 
are computed numerically. The result is 

N N 

^ -^nm) 

n = I m = I 

where the constants ^nm l^e values of the definite integrals of e„, and and their 

derivatives in Eq. (38). The value of I is minimized by picking the values of such that 

N 

a/ ^ 

7 - y [Bnm +A^„)] (Xm-O , 1,2,3, ... ,7V 

da„ ^ 

m = 1 

in matrix form, this equation is 


[75 + B^] [a] - [A +A'^\ [a] = 0 


(40) 


The two solutions of Eq. (40) are ( [a] | = 0, which means the liquid has no motion at all, a solution 
which is not of interest here, and 


I [5 +5^] -a^[A =0 (41) 

Eq. (41) represents a standard eigenvalue problem, and the values of fZ = which allow the determinant 
in Eq. (41) to be equal to zero are the natural frequencies or slosh frequencies, of the problem. Once 
the are known (there are N of them), the constants which are the eigenvectors, can be computed 
with the aid of Eq, (40). 


In the limit as 4) and e, when determined in the above way, satisfy the boundary conditions, 

Eqs. (27), (28), and (30); also <J> satisfies the incompressibility condition = 0 [Eq. (26)] . Thus, is 
the true velocity potential of the low-gravity sloshing problem, and the slosh mode shapes, forces, moments, 
etc. can be calculated. 
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B. Trial Functions and Results 


To make the Rayleigh-Ritz technique feasible, it is usually necessary to pick trial functions that 
“suit” the problem. Two different sets of functions were evaluated during the present work. The first was 
the set of polynomials used by Rattaya [3] : 

n 

= cos 0 ^ (41) 

k = 1,3,... 

where, in order to satisfy = 0, it is necessary that 


,(") = — 


(/I — A: + 2)(n — k + \) 


^k-2 


A: = 3,5, . . .n 


and 


^(«) = \ for all n 

Rattaya reported excellent results using Eq. (41), with < 12, in his analysis of sloshing in ellipsoidal 
tanks, which neglected low-gravity effects. However, these functions gave very poor results during the 
present study and were abandoned rather early.* 

The second set of functions used in the present study were those that describe sloshing in a cylin- 
drical tank: 


= cos e • J, (\„R) • 


cosh [X^(Z+// + L)] \ 

cosh [X^C/^max + ^+^)] J 


(42) 


where =F(R^). 


The constant in each was picked so that satisfied 
in some sense the boundary condition d<t>ldn = 0 on the tank walls; 
this should improve the rate of convergence of l^^e true 

d>. One method that was tried was to use the value of X^ that made 
= 0 dt R = 1. Thus, with this method, the normal velocity 
was zero over all the cylindrical surface, R = i,L— L~ 
//, 0 < 0 < 277, and 0 < ^ < 1, Z = - L - //, 0 < 0 < 277, which 
bounds the actual tank; this is the technique suggested by Moiseev 
and Petrov [8] . (Note, however, that the assumed contact line con- 
dition e = 0, which should be satisfied by the sum is not 

satisfied by each individually when computed in this way; the 
reason for this is that the condition used to calculate X„ does not 
insure that = 0 at the contact point.) Table I gives repre- 

sentative values of for a spherical tank that were calculated using 
these X„, for the case Nqq = 100 and a free surface that is not bent 
over. As can be seen, the comparison with the present results and 
those given by Rattaya [3] (for Nqq - °^) are fairly close; with 
further refinement in the numerical work, it is believed that the 
comparison would be even closer. 


TABLE I. VARIATION OF AS A 
FUNCTION OF TANK FILLING 
AND NUMBER OF TERMS IN 
APPROXIMATING SERIES 
FOR = 100 


Liquid 

Volume 

N 

j 

n 1 , Ref. 3 

47.2% 

5 

12 

1.274 

1.276 

1.230 

35% 

5 

12 

1.239 

1.239 

1.185 

30% 

5 

12 

1.231 

1.231 

1.160 

25% 

5 

12 

1.226 

1.226 

1.140 

20% 

5 

12 

1.225 

1.225 

1.125 

15% 

5 

12 

1.230 

1.230 

1.100 


^Because Rattaya neglected surface tension and curvature effects, he was able to include the boundary condition Eq. (21) 
in his variational integral. Thus, for any finite value ofN^Q, his integral and the present one cannot be made to be identical. 
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The second method was to use the value of in Eq. (42) that made = 0 at the contact line; 

this automatically forced each to be equal to zero at the contact line, which thus identically satisfied 
one more of the boundary conditions. The results using this method were quite similar to the results 
presented in Table I. 

For bent-over free surfaces, neither method of picking X,„ seemed to give satisfactory results. As an 
example, the case Nqq = 100 for a 70% full spherical tank resulted in = 1.797 for \!/ = 0° whereas, 
according to Refs. 5 or 6, the correct result is about = 1.34 when ^ = 5°. After many more test cases 
and manipulations of the trial functions , with no further improvement in the prediction of , it was 

concluded that either (1) a suitable set of functions had not been discovered* or else (2) the boundary 
condition that = 0 at /? = 7?^^ forced the frequency to be too large because of the added “stiffness” 
given to the free surface as a result of the “stuck” contact line. A third alternative is that the present result 
for i// = 0° is correct, or could be made so with further refinement of the technique, so that therefore the 
difference between the present and the result given in Ref. 5 is due to \p being 5° in Ref. 5. This 
alternative would be a likely explanation if the numerical technique used in Refs. 5 or 6 does not converge 
to a definite result as ^ 0°. However, it was not feasible to evaluate the numerical convergence of these 
methods during this project. 


Another obvious set is Legendre polynomials. 



V. CONCLUSIONS 


It has been shown that Hamilton’s principle leads to a variational form of the sloshing equation which 
is equivalent to the usual differential equations. A second form of the variational integral, which may be 
more convenient to use in numerical work, was derived by direct integration of the differential equations 
and boundary conditions and was used in the subsequent work. 

The “contact line condition” implicit in the variational form of the sloshing equations (or assumed in 
the differential equation form) was shown to be ambiguous for the case of a zero-degree contact angle. The 
ambiguity is a result of both the linearization of the equations with respect to the wave height and the 
impossibility of mathematically defining the slosh height at the walls of an ellipsoidal tank. Although the 
ambiguity does not appear to be serious for non-zero contact angles, it is believed to be a critical limitation 
in the analysis of bent-over free surfaces with a zero-degree contact angle. 

Numerical results of sloshing in ellipsoidal tanks for free surfaces that are not bent over (tanks less 
than half-full) compare well with previous work. For bent-over surfaces, however, the computed natural 
frequencies were considerably larger than those computed previously for a contact angle of 5 degrees. The 
discrepancy is probably the result of the contact line condition difficulty mentioned above which, for a 
zero-degree contact angle, implies analytically that the slosh wave is “stuck” to the contact line. A satis- 
factory resolution of the contact line difficulty was not obtained. 
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APPENDIX 


Summary of 

"Low Gravity Sloshing in Rectangular Tanks" 
by 

Franklin T. Dodge 
Luis R. Garza 


Technical Report No. 1 
Contract NAS8-24022 


January 1970 



Liquid sloshing in rectangular tanks was studied theoretically and experimentally under low Bond 
number conditions. The static-free surface shape was computed accurately by an approximate technique, 
and the results were used in the equations of motion for the fluid to determine the sloshing parameters. The 
sloshing equations were solved by Galerkin’s method. The natural frequency parameter was found to 
increase, and the equivalent slosh mass to decrease, under low Bond number conditions. Nonlinearities in 
the experimental results prevented a close comparison of theory and test, but the trends of both were 
similar. Exploratory tests with square tanks showed that nonlinear effects prevailed here also, even for large 
Bond numbers. 

Figures A-1, A-2, and A-3 show the comparison of test and theory, for rectangular tanks having a 
width-to-thickness ratio of 10. 



0 5 10 15 20 25 30 35 40 

BOND NUMBER. NgQ = />gw*/a 


FIGURE A-1. NATURAL FREQUENCY VS FOR DEPTH RATIO OF 1.5 
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BOND NUMBER, N^Q-pg^^lo 


2574 


FIGURE A-2. NATURAL FREQUENCY VS N^q, FOR DEPTH RATIO OF 0.75 



FIGURE A-3. NATURAL FREQUENCY VS FOR DEPTH RATIO OF 0.50 
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